Exact solution for the degenerate ground-state manifold of a strongly interacting 

one-dimensional Bose-Fermi mixture 
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We present the exact solution for the many-body wavefunction of a one-dimensional mixture of 
bosons and spin-polarized fermions with equal masses and infinitely strong repulsive interactions 
under external confinement. Such a model displays a large degeneracy of the ground state. Using a 
generalized Bose-Fermi mapping we find the solution for the whole set of ground-state wavefunctions 
of the degenerate manifold and we characterize them according to group-symmetry considerations. 
We find that the density profile and the momentum distribution depends on the symmetry of the 
solution. By combining the wavefunctions of the degenerate manifold with suitable symmetry and 
guided by the strong-coupling form of the Bethe-Ansatz solution for the homogeneous system we 
propose an analytic expression for the many-body wavefunction of the inhomogeneous system which 
well describes the ground state at finite, large and equal interactions strengths, as validated by 
numerical simulations. 

PACS numbers: 05.30.-d,67.85.-d,67.85.Pq 



I. INTRODUCTION 



Ultracold atomic gases provide a versatile and con- 
trolled system for the study of quantum correlations and 
fluctuations which are particularly strong in one dimen- 
sion (ID). Experiments on two-dimensional optical lat- 
tices [lHj] or on a chip trap [5| have reached the strongly 
interacting Tonks-Girardeau regime. In such impenetra- 
ble boson limit, repulsive interactions play the role of the 
Pauli exclusion principle and the many-body wavefunc- 
tion can be exactly obtained by mapping onto the one of 
noninteracting fermions [6|. Experimental advances on 
trapping and cooling ultracold Bose-Fermi mixtures 0- 
OjJ and the possibility of trapping both species in tight 
atomic waveguides have boosted a theoretical activity 
on ID mixtures. At increasing boson-fermion repulsions, 
mean-field [13j and Luttinger liquid analysis at weak cou- 
pling 14] predict an instability towards phase separation 
(i.e. demixing) of the two components. For a highly sym- 
metric model with equal masses and coupling constants, 
further progress can be made by use of exact solutions. 
For the homogeneous system, a Bethe-Ansatz solution 
is known [15|-[l7| and no demixing is found. The long- 
wavelength properties of its correlation functions have 
been studied using conformal field theory [l8j . Inhomo- 
geneous systems, as in the case of experiments, bring 
about novel issues, such as the spatial structure of the 
ground state. At intermediate interaction strength a 
partial demixing of the two clouds has been found by 
a local density approximation on the Bethe-Ansatz so- 
lution both at zero and finite temperature [l?], EH, H3| • 
In the Tonks-Girardeau limit of infinitely strong boson- 
boson and boson-fermion repulsions a large ground state 



degeneracy is expected [2l| , and is associated to the free- 
dom of fixing the sign of the many-body wavefunction 
under the exchange of a boson with a fermion. For an 
inhomogeneous system, one exact solution of the degen- 
erate manifold has been proposed in 2l| and analyzed 
in detail in [22|-[24j]. The corresponding density profiles 
display no demixing among the two species. Till now no 
expression was known for the other wavefunctions of the 
manifold. In this work we solve several open theoretical 
issues. First of all we find an exact analytical solution for 
all the wavefunctions of the degenerate manifold in the 
Tonks-Girardeau limit, thus generalizing the solution of 
[2ll ] . Secondly, we characterize the solutions in terms 
of their symmetry properties according to group theory 
considerations. At difference from fermionic or bosonic 
spinor systems [25|, [26|, where the state of the system 
can be labelled on the basis of the spin quantum num- 
ber, in order to label the states of the Bose-Fermi mixture 
we introduce a suitable Casimir operator which reflects 
the mixed symmetry under particle exchange. Further- 
more, we find that such symmetry considerations allow 
for the understanding of the shape of the momentum 
distribution, which depends on the choice of the wave- 
function within the manifold. Finally, by linear combi- 
nation of the basis wavefunctions of the degenerate man- 
ifold we individuate the wavefunction which corresponds 
to the ground state at finite, large and equal interactions 
strengths, and we confirm this prediction by comparing 
with numerical DMRG simulations. For the nontrivial 
case of the Bose-Fermi mixture with large degeneracy, 
this analysis allows for the first time to draw a link be- 
tween the Bethe-Ansatz solution of the homogeneous sys- 
tem and the Tonks-Girardeau solution of the inhomoge- 
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neous system. Our solution sheds light onto the general 
properties of the ground state wavefunction of a fully 
quantum problem in the strongly interacting limit. 



II. ORTHONORMAL BASIS SET FOR THE 
DEGENERATE GROUND-STATE MANIFOLD 

A. General considerations 

We consider the model of Nb bosons and Nf spin- 
polarized fermions of masses raj = mp — m, conhncd 
by the same external potential. The particles interact 
via the contact potentials, vbb(x) — gBB$(x), vbf(x) — 
gBFS(x), and we focus on the limit qbb = 9bf — > oo. 
The effect of contact interactions can be replaced by the 
boundary condition that the wavefunction vanishes at 
each BB or BF contact, i.e. 

^(.., Xj, .., X£, ..) =0 whenever Xj — X£. (1) 

We adopt the convention that {xi, xn b } are bosonic 
coordinates and {xn b +x, ...,Xn} are fermionic ones. The 
ground state has a large degeneracy C- B — N\ /Nb ! /Np ! , 
which can be interpreted as choosing Nb positions for the 
bosons out of N = Nb + Np, and amounts to fixing in 
several possible ways the sign of the wavefunction under 
the exchange of bosons with fermions. 

In order to determine an orthonormal basis set for the 
degenerate manifold we proceed as follows. Consider a 
fermionic Slater determinant made of the first total N 
orbitals, 

9 P (xx,...,x N ) = -^=det|&(a;/)], (2) 

where j, £ — 1, ...,iV, and 4>j(x) are obtained by the 
solution of the single-particle Schroedinger equation in 
the potential V(x). ^f displays the correct nodes at 
each BB and BF contact. In a given coordinate sector, 
Xpn) < xp(2) < ■•■ < Xprm, with P being a permutation 
among the N particles, the required many-body wave- 
function is proportional to ^ p(xi, ...,xn)- A useful set 
of orthonormal basis is given by the "snippets" [25[ 

(x u ...,x N \P) = y/Nl\9 F (xi,...,x N )\ (3) 

in the coordinate sector xpm < xpi^) < ■■■ < ^p(iv), 
and zero otherwise. To build the required basis for the 
manifold, we now collect the snippets which correspond 
to exchanging only the positions of the bosons or of 
the fermions among themselves, using the Fermi or Bose 
statistics to fix the relative sign of the various terms. The 
number of groups of snippets subdivided in such a way is 
exactly C- B . This yields the required orthonormal basis 
set {^ a } for the degenerate manifold, since each snippet 
is orthogonal to another and is used only once. 



The density profiles associated to each wavefunction 
ty a are given by 

n B (x) = Nb J dx2---dxN\^a(x,X2--XN)\ 2 , 

n F {x)=N F J dxi...dx N -i\^ a {xi, .., Zjy-i, z)| 2 , (4) 
which is equivalent to computing 

N 

n B(F)( x ) = Y.Pi,B(F)Pi( x )^ (5) 

i=l 

with ViB(F) = 1 if a boson (fermion) is at position i = 

1, N in the configuration a and zero otherwise, and 

Pi(x) = / dx 1 ...dx N \ty F (xi, ..x N )\ 2 6(x - Xi). (6) 

J Xi<X2--.<X N 

B. An illustration with Nb = Nf = 2 

We illustrate the idea in the case N — 4, Nb = Np = 

2, and take a harmonic confinement V(x) = \rwJ 1 x 1 for 
simplicity, assuming the same trapping frequency to for 
the two species. We denote the bosonic coordinates by 
xi, X2 and the fermionic ones by £3, X4, and we expect a 
six-fold degeneracy. We label the basis set using the po- 
sitions of the particles, i.e. BBFF, BFBF, BFFB, FBBF, 
FBFB, FFBB. Let us call this basis the "BBFF" basis. 
According to the above prescription the first wavefunc- 
tion is 

^bbff= \[{xx,x 2 , a; 3 ,a;4|(e+(12))(e- (34)))], (7) 

where by (j£) we denote the permutation between the 
particles j and £, e is the identity permutation, and we 
adopt the usual convention of product among permuta- 
tions [27| ; the snippet basis satisfies (x%, ..Xn\P + Q) = 
(xi, ..xn\P) + (xi, ..xn\Q) . Similarly, the second wave- 
function is obtained by 

Vbfbf = ^[(x 1: x 3 ,x 2 ,x 4 \{e+ {12)){e- (34)))]. (8) 

The other wavefunctions are built in the same way, taking 
as initial coordinate sector the one where the bosonic and 
fermionic coordinates are each in ascending order. 

We display the density profiles of the six basis states in 
FigLT] Each peak corresponds to the position of a particle 
in the BBFF sequence, hence the basis set recalls the 
one of distinguishable particles, as in the case of spinor 
bosons [25(. 

III. SYMMETRY CHARACTERIZATION 

A. Casimir invariance of the manifold 

We would like to label the basis vectors by some ad- 
ditional quantum number. We proceed by exploiting the 
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FIG. 1. (Color online) Bosonic (magenta solid line) and 
fermionic (blue dashed line) density profiles (in units of a^ 1 = 
\J raw / h) as a function of the spatial coordinate x (in units 
of a,ho) for each BBFF wavefunction of the Nb = Nf = 2 
mixture. The total density (black solid thin line) is included 
in the first frame for reference. 



FIG. 2. (Color online) Bosonic (magenta solid line) and 
fermionic (blue dashed line) density profiles (in units of a^ 1 ) 
as a function of the spatial coordinate x (in units of aho) with 
the Y (top) and Y' (bottom) symmetry for TVs = Nf = 2. 



exchange symmetry between bosons or fermions among 
themselves. According to a group theoretical analysis, 
there are only two possible Young tableaus associated 
to a quantum mechanical system of mixed Bose-Fermi 
symmetry, i.e. symmetric in its first Nb coordinates and 
antisymmetric in its last Np: 



Y 



Fi 






Bn b 






Fn f 



Y' 



Bx 




Bn b 


Fi 






Fn f 



, (9) 



each with a dimension of C 



C 



N-l 

N F 



respectively. 



To each Young tableau it is possible to associate a value 
of a Casimir invariant, obtained from the generators of 
the permutation group Sn, and which commutes with all 
elements of the group. In this particular case we choose 
C = J2i<j(ij)i the sum over au tne transpositions [28| . 
Two different eigenvalues of the Casimir operator are as- 
sociated to the two tableaus, namely 



Cy 



N B + 1 
2 

N B -1 



N F 



N F + 1 



(10) 



with degeneracy corresponding to the dimension of the 
tableau. We remark that the eigenvalue Cyry) essentially 
counts the number of two-particle exchange allowed by 
the corresponding tableau, where a symmetric permuta- 
tion (from the row) is counted as +1 and an antisymmet- 
ric permutation (from the column) is counted as —1. 

We represent the Casimir operator on the BBFF or- 
thonormal basis set (see appendix |A"1 for some examples). 
The eigenvectors of the Casimir operators, obtained by 
diagonalization, are expressed as linear combinations of 
the BBFF basis vectors, and provide an alternative (but 
not orthonormal) basis set ^(sci, xjy) characterized 
by the symmetry of the tableau. 



The related bosonic and fermionic momentum distri- 
butions 



i(F)(P) = ^j dx dx'e-^-^p^ B(F) (x,x') (11) 



can be computed from the bosonic and fermionic one- 
body density matrix 



B. 



p B (x,x') = N B J dx2---dx N ^^(x,x 2 , ...,xn) 

X^!j,(x',X 2 , ...,X N ), 

p F (x,x') = N F jdxi...dx N -i^(xx, ,.x N -i,x) 

x* m (xi,...xat_i,x'). (12) 



Density profiles and momentum distributions at 
a given symmetry 



We return to the earlier example of two bosons and two 
fermions in a harmonic confinement. The density profiles 
of the eigenstates of Y and Y 1 are shown in Fig. [5J We 
notice that in this highly symmetric problem the role of 
bosons and fermions is simply exchanged among the two 
submanifolds. 

We also display the corresponding momentum distri- 
butions in Fig J3J It is worth mentioning that the momen- 
tum distributions associated to the Y symmetry display 
less peaks (with a number of peaks in np{p) equal to Np 
as in [22j]) and are narrower than those associated to the 
Y' symmetry, following the intuition that the symmetry 
of the latter tableau is more "Fermi-like" (i.e. extended 
in the vertical direction). 
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FIG. 3. (Color online) Bosonic (magenta solid line) and 
fermionic (blue dashed line) momentum distributions (in units 
of p^o — l/V hmijj) as a function of p (in units of pho) with 
the Y (top) and Y' (bottom) symmetry for Nb = Nf = 2. 
The first and second panels in each row coincide. 



IV. GROUND STATE AT LARGE, FINITE 
INTERACTION STRENGTH 

A. Analysis of special solutions 

At finite interactions gss = gBF, the ground state 
is expected to display the Y symmetry [l5|, [29|, since 
the associated wavefunction has less nodes than the one 
with Y' symmetry. From a continuity argument starting 
from the noninteracting solution, we also expect that it 
is nonvanishing in all permutation sectors. In Ref.f2lj a 
special solution with the latter property was proposed, 

*gm = A bb A bf ^f{xi,... ,xn), (13) 

where the mapping functions are 

Abb = Hi<j<e<N B sgn(xj - x t ), 

A B f = Hi<j<N B ,N B +i<e<N sgn(xj - x t ). (14) 

We test its symmetry by evaluating the average value 
of the Casimir operator (C) = (*|C|*)/(*|*). The ex- 
pression of the (unnormalized) GM wavefunction on the 

BBFF basis * GM = Y^f pi/NBl c ^ a is c a = 1 for 
any a. From the representation of the Casimir operator 
(see appendix [AJ we obtain that for N B = Np = 3 the 
average value of the Casimir operator corresponds to its 
maximal eigenvalue Cy — 3, hence the GM wavefunction 
has the symmetry of the Y tableau. This is not the case 
for N B — Np — 2 where a wavefunction which is nonva- 
nishing in all coordinate sectors and has the Y symmetry 
is ^' GM = A BB ^p(xx, ...,xn). More generally, we have 
proven that the GM wavefunction with odd Nb = Np 
has always the symmetry of the Y tableau (see appendix 
|B1 for demonstration) . 

For the homogeneous system, the Bethe-Ansatz 
method provides a solution of the model of a Bose-Fermi 
mixture with equal bosonic and fermionic masses and fi- 
nite but equal coupling strengths qbb = 9bf [SIB El- 
The solution for the many-body wavefunction is built 



with the symmetry of a given tableau [3fj, l31| . In anal- 
ogy with fermionic systems, the solution is expressed in 
terms of spatial coordinates and "pseudospin" integer co- 
ordinates yi which correspond to the relative positions of 
the bosons in the coordinate sector Xpt\\ < Xp(2) < ... < 

xp(m, namely {yi,...y NB } = {P^i 1 ), -P -1 ^)}- In 
Ref.[19] Imambekov and Demler study the strongly inter- 
acting limit of the Bethe-Ansatz solution (BA) for N B , 
Np odd. They notice that the wavefunction decouples as 
a product of an "orbital" part and a "spin" part, as 

^>ba ~ det[e i % K ^}* F (x 1 ,...x N ), (15) 

where the orbital part is a Slater determinant of the first 
N orbitals, of course chosen for the homogeneous prob- 
lem V F (xi,...x N ) = (l/VM)det[e^^], and the "spin" 
part includes only the bosonic coordinates j/i, ...t/iv B , 
with k = {-(N B - l)/2 + N/2, N/2, ...(N B - l)/2 + 
N/2} for the ground state. Taking advantage of the 
decoupling among orbital and spin part of the wave- 
function, we generalize such a solution to the inhomo- 
geneous system, by replacing the orbital part of (|T5|) by 
its corresponding expression under confinement, Eq. ([2]l, 
We check its symmetry by evaluating the average of the 
Casimir operator. For Nb — Np = 3 an explicit calcula- 
tion by expansion of (fT5j) on the BBFF basis yields the 
maximal value (^>ba\C\^>ba)/(^ba\^ba) = 3, implying 
that the generalized BA wavefunction has the Y symme- 
try. By construction, the BA solution for the trapped 
case has the same form as the one obtained for arbitrary 
interactions in the homogeneous system [34[ . Notice that 
although the GM and generalized BA wavefunctions are 
not simply proportional to each other, their density pro- 
files coincide, displaying no demixing. We remark that 
Eq. (|15|) can generate a wavefunction with the Y' sym- 
metry by adopting a different choice of spin rapidities, 
e.g. for K = {-(AT B -l)/2,...,0,...(AT s -l)/2} we obtain 
an average Casimir operator equal to —3. 

It is important to notice that the form of the solution 
for the wavefunction at finite large interactions depends 
from the way in which the limit 555—7-00 and gBF ^ 00 
is approached. Consider for example the case where gsB 
tends to infinity and gsF is finite. The above symmetry 
classifications are not useful in this case. The bosonic 
component can be mapped onto a fermionic one, and the 
problem can be reduced to the one of spin 1/2 fermions 
as in 1 35] . Also in such decoupling of spatial and 

spin degrees of freedom is predicted with a different form 
for the spin part with respect to Eq. (fT5| . 



B. Numerical illustration 

We have tested our predictions by comparing with nu- 
merical DMRG simulations (3(|. The DMRG techniques 
provides an efficient numerical solution for the lattice 
model of an interacting Bose-Fermi mixture subjected 



5 




FIG. 4. (Color online) Bosonic (magenta solid line) and 
fermionic (blue dashed line) density profiles as obtained from 
DMRG simulations (in units of a^ 1 ) as a function of x 
(in units of a^o) at increasing interaction strength U/t = 
1,10, 10 2 ,10 4 (from left to right and from top to bottom). 
The harmonic trap strength is V/t = 7 x 10~ 6 and the num- 
ber of lattice sites used in the simulation is L = 128. In the 
last panel the analytical prediction from the generalized BA 
wavefunction (thin black solid line) is shown. 

to an external potential, described by the Hamiltonian 

3 3 

+U B b n B ,j(n B ,j - 1) + U B f n B , n F ^ 

3 3 

+ J2V(J - L/2) 2 (n Bj +n Fj ).(16) 

3 

Here, bj and fj are the bosonic and fermionic field op- 
erators acting on a site j, with corresponding density 
operators nsj — bjbj and riF.j — flfj, t is the tunnel 
constant, Ubb and Ubf are the on-site BB and BF in- 
teractions, V is the strength of the external harmonic 
confinement, which is taken to be the same for bosons 
and fermions, and L is the number of sites in the lattice. 

Since in the limit of very low lattice filling the lattice 
model reproduces the continuum model, we have used 
this method to obtain the ground-state density profiles 
and momentum distributions for the Bose-Fermi mixture 
at increasing interaction strength. For the sake of com- 
paring with the previously proposed special solutions we 
have restricted ourselves to the case of equal BB and BF 
interactions, U = Ubb = Ubf- Figure 0] displays the 
bosonic and fermionic density profiles for a Bose-Fermi 
mixture with Nb — 2, Np — 2 at increasing interaction 
strength. The density profiles evolves from a partially 
demixed one at intermediate interaction strength, to a 
nondemixed one at large interactions, in agreement with 
the predictions of the GM and generalized BA wavefunc- 
tions. This behaviour was also noticed in a density func- 
tional calculation [32| . For a Bose-Bose mixture a similar 
absence of demixing at strong coupling was observed [33j . 




FIG. 5. (Color online) Bosonic (magenta solid lines) and 
fermionic (blue dashed lines) momentum distributions (in 
units of pj^ — l/^hmu) as a function of p (in units of 
Pho), from the numerical DMRG data at increasing interac- 
tion strength (U/t = 1, 10, 10 2 , 10 4 , from top to bottom and 
from left to right, thick lines). At increasing interactions the 
curves approach the analytical predictions from the general- 
ized BA wavefunction for the bosonic (thin black solid lines) 
and fermionic (thin black dashed lines) momentum distribu- 
tions, shown in each panel for reference (almost indistinguish- 
able from the numerical data in the last panel). The other 
parameters used in the simulation are V/t = 7 x 10 -6 and 
L = 128. 



The momentum distributions for the bosonic and 
fermionic components are illustrated in FigjS] as com- 
pared with the predictions of the generalized BA wave- 
function. At increasing the interaction strength the mo- 
mentum distributions approach those obtained by the 
generalized BA wavefunction, showing that this solution 
accurately describes the trapped mixture at finite, large 
but equal interactions. 



V. OUTLOOK AND PERSPECTIVES 

In this work we have studied the exact solutions of a 
highly symmetric model for Bose-Fermi mixture in the 
strongly interacting limit. Our results are relevant for 
the ongoing experiment on ultracold mixtures of atomic 
gases in tight atomic waveguides, with particular atten- 
tion of the case of 173 Yb- 174 Yb Bose-Fermi mixtures [l2| 
where the fractional mass difference among the two iso- 
topes is small. We have provided first an orthogonal basis 
set of solutions which span the degenerate manifold. Sec- 
ondly, we have grouped such solutions on the basis of the 
Casimir invariant, associated to a given Young tableau. 
Finally, we have analyzed two special solutions of the 
problem and discussed their symmetry according to the 
average value of the Casimir operator. By comparing 
with DMRG simulations, we have found that the wave- 
function obtained by generalizing the Bethe-Ansatz solu- 
tion to trapped systems accurately describes the density 
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profile and momentum distribution of the mixture at fi- 
nite, large but equal BB and BF coupling strengths. This 
wavefunction can be used to describe the Bose-Fermi 
mixture in arbitrary external potential. Splittings and 
mixing of the states with different symmetries are ex- 
pected when different masses and different BB and BF 
coupling constants are chosen. Our solution serves as a 
guideline for further numerical studies. This work opens 
also the way to the study of the dynamical properties of 
the strongly interacting Bose-Fermi mixture. Signatures 
of strong correlations could be found in the collective ex- 
citation spectrum. It would also be interesting to inves- 
tigate how particular states in the degenerate manifold 
can be prepared and addressed. 
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Appendix A: Casimir operator on the BBFF basis 

For Nb = 2, Np = 2 the representation of the Casimir 
operator C = X)i<j(i?) on the orthonormal BBFF con- 
figuration basis reads 

/ 1 -1 1 -1 \ 
1 110-1 

- 1 1 1 - 1 . (Al) 
110 11 

-10 110 1 

\ -1 -1 1 1 / 

For Nb = 3, Np = 3 the Casimir operator reads 
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(A2) 



Appendix B: Demonstration of the symmetry of a 
special wavefunction 

We prove here below that for odd Nb = Nf the 
GM wavefunction has the symmetry of the Y tableau. 



This readily follows from the simple expression of 
the GM wavefunction on the BBFF basis ^gm — 
^2N\/N F \/N B \ an( j f rom t ne f ac t that for odd values of 
Np and Nb each line of the matrix which represents the 
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Casimir operator on the BBFF basis has N B Np non- 
vanishing entries with value 1 in N B (Np + l)/2 cases 
and value -1 in N B (Np — l)/2 cases. Each line of the 
matrix corresponds then equally to the value of the av- 
erage of the Casimir operator, (C) = (*|C|*)/(*|*), 



and its value corresponds to the maximal eigenvalue 
Cy = N B = N F . For other values of N B , N F the repre- 
sentation of the Casimir operator is more complicated to 
predict and the above demonstration does not hold. 
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